Pseudo Parity-Time Symmetry in Optical Systems 
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We introduce a novel concept of the pseudo parity-time iVT) symmetry in periodically modulated 
optical systems with balanced gain and loss. We demonstrate that whether the original system 
is PT-symmetric or not, we can manipulate the property of the VT symmetry by applying a 
periodic modulation in such a way that the effective system derived by the high-frequency Floquet 
method is VT symmetric. If the original system is non-PT symmetric, the VT symmetry in the 
effective system will lead to quasi-stationary propagation that can be associated with the pseudo 
VT symmetry. Our results provide a promising approach for manipulating the VT symmetry of 
realistic systems. 
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Parity-time {VT) symmetry is an important concept 
in physics recently developed in application to optical 
systems. In quantum mechanics, the fundamental re- 
quirement of Hermitian Hamiltonians guarantees the ex- 
istence of real eigenvalues and probability conservation. 
However, it was found recently that a wide class of non- 
Hermitian Hamiltonians with the (parity-time) VT sym- 
metry can still possess entirely real eigenvalue spectra ,1- 
Although the extension of quantum mechanics based 
upon non- Hermitian but 7'7~-symmetric operators is still 
a subject to debates, optical systems with complex refrac- 
tive indices 5-d2j are widely used to test the VT sym- 
metry in non-Hermitian systems, because of the equiva- 
lence between the Schrodinger equation in quantum me- 
chanics and the wave equation in optics [l3[. In the 
last few years, the VT symmetry was observed in sev- 
eral optical systems with gain and loss, such as optical 
3, [l^ , microwave billiard pJi] , and large-scale 



.5^ . couplers [14 
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temporal lattices [17 1 



Similar to the electron transport in periodic crystalline 
potentials and the quantum tunneling in periodically 
driven systems [3] , the light propagation in waveguides 
can be effectively controlled by periodic modulations [isil - 
[23} . In an optical system, periodic modulation is as- 
sociated with a periodic refractive index. Mathemat- 
ically, an optical system with periodic complex refrac- 
tive index is equivalent to a time-periodic non-Hermitian 
quantum system. It has been theoretically demonstrated 
that high-frequency diagonal modulations can effectively 
rescale the natural coupling strength 



the periodically modulated system may be described by 
a 7^7~-symmetric Hamiltonian [23| . Naturally, an impor- 
tant question arises: Can the VT symmetry appears in an 
effective system even if the periodically modulated system 
is non-VT- symmetric? In other words, can we employ 
periodic modulations to manipulate the VT symmetry? 

In this Letter, we study the light propagation in a pe- 
riodically modulated optical coupler with balanced gain 
and loss and apply a bi-harmonic modulation along the 
propagation direction. The Hamiltonian for the modu- 
lated system is non-PT-symmetric if the relative phase 
between the two applied harmonics is not or n. Ap- 
plying the high-frequency Floquet approach, the modu- 
lated system is effectively described by an effective aver- 
aged system, whose VT symmetry can be manipulated 
by tuning the modulation amplitude/frequency. More 
importantly, the VT symmetry can appear in the effec- 
tive system corresponding to a non-'PT-symmetric and 
non-Hermitian Hamiltonian. Different from the VT sym- 
metry from a "PT-symmetric Hamiltonian, which leads 
to stationary light propagation of bounded intensity os- 
cillation, the VT symmetry from a non-PT-symmetric 
Hamiltonian will lead to quasi-stationary light propaga- 
tion of unbounded intensity oscillation. Therefore, we 
term the induced symmetry associated with modulated 
systems as the pseudo VT symmetry. 

In optics, the electric field E{x, z) of light obeys the 
wave equation. 
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periodically modulated Hamiltonian can be mapped into 
a time-independent one. Like the unmodulated system, 
the VT symmetry may appear in the effective system if 



.dE{x, z) 
dz 



1 d^E{x, 
2k dxF 



+ V{x,z)E{x,z), (1) 



where, k = k^riQ, fco = 27r/A and V{x, z) — fco["-o ~ n{x)\ 
with the substrate index no, the free-space wavelength 
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A and the complex refractive index distribution n{x) = 
TT-o + nii{x,z) + ini{x), where n/j and n/ are real and 
imaginary parts of the refractive index. Therefore, the 
effective potential reads as V{x, z) — VR(a;, z) + iVi{x) = 
— fco[^i?(2;, z)+ini(x)]. With the ex per imental techniques 
developed in recent years (gI. [lol. [l2l.1l^. [isj , one can make 
Vi{-x) = -Vi{x) and Vr{x,z) = Vb(x) + Vi(x, z) with 
the unmodulated part Vq{x) being a symmetric double- 
well function and the modulation Vi{x,z) ~ V'{x)F{z) 
described by an anti-symmetric function V'{—x) — 
—V'{x) and a bi-harmonic function see Fig. 1. 

Using the coupled-mode theory, the electric field for 
a two-channel coupler can be expressed as E{x, z) ~ 
ci{z)ipi(x) + C2{z)'ip2ix) with the two transverse dis- 
tributions {^pi{x),'ip2ix)} and the two field amplitudes 
{ci{z),C2{z)}. Integrating Eq. (1) with respect to x, we 
find {ci(z), 02(2:)} obeying 



transformation 
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, (2) 



with the dimensionless propagation distance z, the inter- 
channel coupling strength w, the gain/loss strength 7 and 
the bi-harmonic modulation 



S{z) = -A[sm{ujz) + f sin(2a;z -f 0)] . 



(3) 



Here, cj) G [0, 2tt) denotes the relative phase between the 
two harmonics, uj is the modulation frequency, A is the 
modulation amplitude and / is a dimensionless coeffi- 
cient. Since the system is invariant under the transfor- 
mation C2 — C2 and v —v, below we will only con- 
sider the case of u > 0. Defining the parity operator 
as P, which interchanges the two channels labeled by 1 
and 2, and the time operator as T: i — —i, z —z, 
which reverses the propagation direction, the Hamilto- 
nian H{z) for the modulated system (2) is VT symmet- 
ric if PTH = HPT. If = or tt, S{-z) = -S{z), 
H{z) is PT-symmetric. Otherwise, if 7^ and tt, 
S{zq — z) ^ —S{zq + z) for any arbitrary constant zq, 
H{z) becomes non-'PT-symmetric. 




FIG. 1: Schematic diagram of a periodically modulated two- 
channel optical coupler with balanced gain and loss. The pe- 
riodic change of color along z-axis denotes the periodic mod- 
ulation. 

Under the condition of w ^ {'y, 7}, one can implement 
the high-frequency Floquet treatment. Introducing the 
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ciexp|— i — cos(cl;z) -I- cos(2wz -I- (/)) |, (4) 
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C2 exp • 



r r A Af 1} 

l+i — cos(ajz) + ^— cos(2wz + 0) |, (5) 



and averaging the high frequency terms, one can obtain 
an effectively unmodulated system 
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with the effective coupling strength 
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exp(im(/)). (7) 



The modulus of J depends on the values of A/oj and <j). 
If A/uj is relatively small, the modulus \J\ is almost in- 
dependent on the relative phase 4>. When A/oj increases, 
the modulus | J| becomes sensitively dependent on (j). In 
particular, the modulus |J| equals to zero at some spe- 
cific values of A/uj (such as A/uj ~ 2.4 and 5.52) and 
(j> = 7r/2 or 37r/2. In Fig. 2 (a), choosing / = 1/4, we 
show the contour plot of | J| as a function oi A/uj and (j>. 

By diagonalizing the Hamiltonian for the effective sys- 
tem (6), the two eigenvalues are given as 



±|J| 



1 - 



7 
2IJI 
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(8) 



Obviously, dependent on the values of 2yj| , the two eigen- 
values can be real or complex. The two eigenvalues are 
real if 7 < 2\J\ and they become complex if 7 > 2| J|. 
Therefore, ^critical = 2|J| is the critical point for the 
phase transition between real and complex spectra in the 
effective system, which corresponds to the original sys- 
tem (2) under high-frequency modulations. The spon- 
taneous T'T'-symmetry-breaking transition takes place 
in the effective model (6) when the imaginary part of 
e changes from zero to nonzero. Surprisingly, unlike 
our conventional understanding, we find that the quasi- 
energies can be real even if the modulated system (2) is 
non-P7"-symmetric (i.e. (f> ^ and tt). 

The parametric dependence of |Im(e)| is shown in 
Fig. 2 (b-e). In Fig. 2 (b-c), we show |Im(£)| as a func- 
tion of 7 and (p for / — 1/4. For small A/uj, such as 
A/uj = 1 in Fig. 2 (b), |Im(e)| is almost independent on 
the phase cj) and the transition from a completely real 
quasi-energy spectrum (|Im(e)| = 0) to a complex spec- 
trum (|Im(£)| ^ 0) take places when 7 increases. Near 
a minimum of | J|, such as A/ui — 2.4, |Im(e)| strongly 
depends on the phase 0, see Fig. 2 (b). In Fig. 2 (d-e), 
we show |Im(£)| as a function of 7 and A/uj for (d) cj) ~ 
and (e) = 7r/2. Near the minima of \J\, such as A/uj ~ 
2.4,5.52,..., |Im(£)| shows significant difference between 



3 



the two cases of </) = and </) = 7r/2. In particular, at the 
minimum points, \J\ vanished to zero for = 7r/2 and 
the corresponding critical value ^critical = 2| J| is reduced 
to zero. Similar to a non-Hermitian system with no mod- 
ulations, the spontaneous 'PT'-symmetry-breaking tran- 
sition (|Im(£)| = ^ |Ini(e)| ^ 0) can be observed by 
tuning the gain/loss strength 7. More interestingly, for 
our modulated system (2) of fixed 7, it is possible to ob- 
serve the spontaneous 'PT'-symmetry-breaking transition 
by tuning the modulation parameters (j) and A/uj, see 
Fig. 2 (c-e). 
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FIG. 2: (color online) The parametric dependence of the ef- 
fective coupling I J| and the imaginary parts of quasi-energies 
|Im(e)|. Top row [(a)]: |J| versus A/lo and </> for / = 1/4. 
Middle row [(b) and (c)]: |Im(e)| versus ^/vr and 7 for (b) 
A/oj = 1 and (c) A/w = 2.4. Bottom row [(d) and (e)]: 
|Im(e)| versus A/u and 7 for (d) = and (e) (f) = n. The 
other parameters for |Im(£)| are chosen as v — 1, uj — 10 and 
/ = 1/4. The white curves are the boundary {^critical = 2| J|) 
between |Im(e)| — and |Im(£)| 7^ 0. 

Based upon the high-frequency Floquet treatment, it 
seems that, whether the modulated system (2) obeys 
a ■pT-symmetric Hamiltonian or not, completely real 
quasi-energy spectrum always appear if 7 < 2| J|. This 
is obviously inconsistent with the previous theory [ij-Q 
which tells us that only PT-symmetric Hamiltonian sys- 
tems can support completely real spectra. So, what really 
happens in the modulated non-Hermitian and non- "PT- 
symmetric Hamiltonian system? 

In general, for systems under temporally periodic mod- 
ulations of arbitrary modulation frequencies, one can use 
the Floquet theory to calculate their quasi-energies. Sim- 
ilar to the Bloch states for systems of spatially periodic 



potentials, the modulated system ([2]) has Floquet states, 
{ci(z),C2(z)} = e~*^''{ci(z), C2(z)}. Here, the propaga- 
tion constant e is called as the quasi-energy, and the two 
complex amplitudes ci(z) and £2(2) are periodic with the 
modulation period T — 27r/w. 

To check the validity of the effective model (6), we 
calculate numerically the quasi-energies for the origi- 
nal model ([2) and then compare the numerical results 
with the analytical formula (8). Under the condition of 
u! — 10, V — 1 and / = 1/4, the analytical formula (8) 
is well consistent with the numerical results for 7 up to 
3. The analytical and numerical values for the quasi- 
energies e are in good agreement and only show tiny dif- 
ference dependent upon the phase (p. As two examples, 
we show Im(e) (the imaginary part of quasi-energy) ver- 
sus 7 for = and = 7r/2 in Fig.|3](a) and (b), respec- 
tively. It clearly shows that the analytical results (red 
lines) agree well with the numerical results (black lines) 
obtained from the original model 1^. Below the criti- 
cal point ^critical — 2| Jj, for systems described by PT- 
symmetric Hamiltonians (i.e. </> = or tt), the numerical 
results confirm the entirely real quasi-energy spectrum, 
see Fig. 3 (c). However, for systems described by non- 
PT-symmetric Hamiltonians (i.e. 0^0 and tt), the nu- 
merical results show the quasi-energies e still have small 
non-zero imaginary parts even below the critical point 
Icriticai = 2|J|, scc Fig. 3 (d). This means that, if the 
original system (2) obeys a non-PT-symmetric Hamilto- 
nian, the entirely real quasi-energy spectrum for the effec- 
tive model (6) does not correspond to perfectly entirely 
real quasi-energy spectrum for the original system (2). 
Therefore, such a PT symmetry in the effective model (6) 
corresponds to a kind of pseudo PT symmetry in the 
original model (2). 

Through numerical integration, we analyze the light 
propagations in the modulated system (2) of differ- 
ent quasi-energies corresponding to different symme- 
tries. The light propagation sensitively depends upon the 
quasi-energies. Stationary light propagations of bounded 
intensity oscillations appear if all quasi-energies are real. 
Non-stationary light propagations of unbounded inten- 
sity oscillations appear if at least one of quasi-energies 
is complex, in which quasi-stationary light propagations 
of slowly varying time-averaged intensities appear if the 
two quasi-energies for the effective system (6) are real. In 
Fig. 4, for u = 1, A = 10, / = 1/4, = 10 and 7 = 0.1 
(which is below the critical value ^ critical) i we show the 
intensity evolution from the initial state of ci (0) = 1 and 
C2(0) — 0. In which, the two intensities Ij{z) = \cj{z)\ , 
the total intensity It{z) — Ii{z) + ^2(2) and the time- 
averaged total intensity If^{z) = ^ ° It{z)dz with 
Ts = 27r/ |Re(e2) — Re(ei)| and Re(ej) being the real 
part of Ej. In short-distance propagations, /i,2(-z) and 
It{z) oscillate periodically and it is hard to see the dif- 
ference between the cases of = and tt/2, see Fig. 
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FIG. 3: (color online) Comparison between numerical and 
analytical results of Im(e), the imaginary part of quasi-energy. 
Upper row: Im(e) versus 7 for (a) <^ = and (b) <j!> = 7r/2. 
Solid lines are for numerical results obtained from the original 
model ((2| and red dashed lines are analytical results given 
by the formula (8) for the effective model (6). Lower row: 
the enlarged regions of (a) and (b) near the bifurcation point 
given by the analytical formula (8). The other parameters are 
= 1, A = 10, / = 1/4 and = 10. 



4 (c-d). However, significant difference appears in long- 
distance propagations. For a PT-symmetric Hamiltonian 
system (2) of = 0, the time-averaged total intensity 
It'"{z) keeps unchanged, see Fig. 4 (a). For a non-T'T- 
symmetric Hamiltonian system (2) oi (f> — 7r/2, the time- 
averaged total intensity slowly increases, see Fig. 
4 (b). The quasi-stationary light propagations of slowly 
varying I^'"{z) is a direct signature of the pseudo VT 
symmetry in the modulated system (2). 



2000 




FIG. 4: (color online) Intensity evolution from the initial state 
of ci(0) = 1 and C2(0) = 0. Upper row : long-distance time- 
averaged intensity evolution for (a) <^ = and (b) <j) = 7r/2. 
Lower row: short-distance intensity evolution for (c) = 
and (d) <f) ~ ("i-- The other parameters are chosen as « = 1, 
A = 10, / = 1/4, a; = 10 and 7 = 0.1. 

Now, we discuss the experimental possibility to ob- 
serve our theoretical predictions. Recently, several VT- 
symmetric optical systems have been realized experi- 
mentally [i, O, E til [il- 

Complex refractive index 
of gain/loss effects can be obtained from quantum-well 



lasers or photorefractive structures through two-wave 
mixing |25| . Periodic modulations can be introduced by 
out-of-phase harmonic modulations of the real refractive 
index 13, 21, 2^ or p eriodic curvatures along the prop- 



agation direction 



or p eriodic curv 

[H [ii m [ii] 



For a short optical 



coupler under periodic modulations, spontaneous VT- 
symmetry-breaking transitions can be observed, whether 
the system Hamiltonian is "PT-symmetric or not. In such 
a system, the light propagation will be periodic and sta- 
ble if 7 < 2| J|, and an instability will be observed if 
7 > 2| J|. The critical point ^critical = 2| J| can be ad- 
justed by controlling the modulation parameters Ajuj^ f 
and 4> in addition to controlling the gain/loss strength 7. 
However, for a long optical coupler under periodic mod- 
ulations, the light propagation below the critical point 
(i.e. 7 < 2|J|) depends on the Hamiltonian symmetry. 
If the Hamiltonian is PT-symmetric, the light propaga- 
tion is periodic and stable, in which the time-averaged 
total intensity keeps conserved. Otherwise, if the Hamil- 
tonian is non-T'T-symmetric, the light propagation is 
quasi-stationary, in which the time-averaged total inten- 
sity slowly changes. 

In summary, we have studied the non-Hermitian 
Hamiltonian systems under periodic modulations and in- 
troduce the concept of the pseudo VT symmetry. If the 
modulated system obeys a VT symmetric Hamiltonian, 
there exists a truly spontaneous PT-symmetry-breaking 
phase transition from a real quasi-energy spectrum to a 
complex one. If the modulated system obeys a non-'PT'- 
symmetric Hamiltonian, although there exists a spon- 
taneous PT-symmetry-breaking phase transition in the 
effective system derived from the high-frequency Floquet 
treatment, there is no truly spontaneous PT-symmetry- 
breaking phase transition in the original system. Cor- 
responding to the real spectrum for the effective sys- 
tem, the original system has a quasi-real spectrum of 
small imaginary parts, which leads to a quasi-stationary 
light propagation of slowly varying time-averaged to- 
tal intensity. This is the pseudo VT symmetry in the 
non-Hermitian system described by a non-'PT'-symmetric 
Hamiltonian. 

In addition to the discovery of the pseudo VT symme- 
try, we believe that our work brings three key advances 
to related fields. Firstly, although the high-frequency 
Floquet treatment can capture most key features, some 
important information (such as the so-called pseudo VT 
symmetry) may be lost. Secondly, periodic modulations 
provide a new route to the observation of the spontaneous 
PT-symmetry-breaking transition. Thirdly, as the inter- 
channel coupling can be effectively switched off by con- 
trolling the modulation and the corresponding intensity 
grows exponentially even for arbitrarily weak gain/loss, 
this exponential growth offers an efficient way to beam 
amplification in optical waveguides. 
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